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ABSTRACT 

We use gravitational microlensing to determine the size of the X-ray and 
optical emission regions of the quadruple lens system Q 2237+0305. The optical 
half-light radius, log(i?i/2,y/cm) = 16.41±0.18 (at Xrest = 2018A), is significantly 
larger than the observed soft, log{Ri/2,soft/cm) = 15.76lo;34 (1.1 — 3.5 keV in the 
rest frame), and hard, log(i?i/2,/iard/cni) = 15.46to;29 (3.5 — 21.5 keV in the rest 
frame), band X-ray emission. There is a weak evidence that the hard component 
is more compact than the soft, with \og{Ri/2,soft/ Ri/2,hard) — O.SO^q^I. This 
wavelength-dependent structure agrees with recent results found in other lens 
systems using microlensing techniques, and favors geometries in which the corona 
is concentrated near the inner edge of the accretion disk. While the available 
measurements are limited, the size of the X-ray emission region appears to be 
roughly proportional to the mass of the central black hole. 

Subject headings: accretion, accretion disks — gravitational lensing: micro — 
quasars: general 



1. Introduction 



AGNs play a crucial role in the evolution of the Universe, but their internal structure 
remains poorly understood. The energy is generated by accretion onto a supermassive black 
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hole (SMBH ) through a disk, but the detailed structur e of the disk (e.g., |K 



aes et ah 2004 



BlaesI 120071 ) , the origin of the non -thermal emissions (iReynolds fc NowakI l2003l ) , and the 
source of any outflows (lElvisI l2000l ) are all topics of debate. Although simple thin accretion 
disk models d escribe some of the observed behavior, they do not reproduce all the observed 
features (e.g., iBlaed 120071 . and references therein). Testing accretion disk models is chal- 
lenging in large part becau se of our inabil ity to resolve the emission regions. Reverberation 
mapping techniques (e.g., iPetersonI 119931 ) strongly constrain the structure of the distant 
emission line regions, but not the continuu r n emission regions beyond the limits a l ready set 



by time variability (e.g., Kelly et al.l 120091 : iKozlowski et al.ll2010l : iMacLeod et al.ll2010l for 



modern models of optical variability). Fortunately, gravitationally lensed quasars provide us 
with a unique tool to zoom in on the structure of AGN and explore their physics in more 
detail - a "telescope" that only becomes more powerful as the target source becomes smaller. 

Microlensing magnification is caused by stars and white dwarfs in the lens galaxy close 
to the line of sight towards the lensed quasar images. It leads to uncorrelated flux vari- 
ations due to the relative motions of the quasar, the lens, its stars, and the observer 
(see the review by IWambsganssI |2006| ). The microlensing signal depends on the struc- 
tural and dynamical properties of the source and the lens. Since the magnification de- 
pends upon the size of the source, the simplest application of quasar microlensing is to 
measure the size of quasar accretion disks. In essence, the amplitude of the microlens- 
ing variability encodes the disk size, with larger disks showing lower variability ampli- 
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Microlensing has also been used to explore the spatial structure of the broad line re- 



will be small due to its large size (see 



gion (BLR) of quasa rs fe.g.. iLewis fc Bellelll998l : lAbaias et al.ll2002l : ISluse et al.l 120071 12011 



O'Dowd et al.l 120101 : iGuerras et al.l |2012[). although the amplitudes of BLR lensing effects 



Mosquera fc Kochanekll2011bl ) 



The largest microl e nsing v ariability amplitu des a re observed in the X -rays, as empha- 
sized by IChartas et al.l ( 12002| ) , iDai et al.l ( l2003l ) , and iPooley et al.l ( 120071 ) , indicating that 
this emission must arise close to the inner edge of the disk. Quantitative microlensing X-ray 



size estimates are now availab le for several systems ( iDai et al.l l2010l : iMorgan et al.l 12012 



Blackburne et al.ll2011bl |2012|), although the sparse data means that the measurements are 
largely upper limits because the likelihood distributions for the lower limits are prior depen- 
dent. 
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The well-known quadruple lensed quasar Q 2237+0305 ( Huchra et al.lll985l ) remains the 
best system for studying microlensing phenomena. In this system a zs = 1.69 background 
quasar is lensed by a nearby spiral galaxy at zl = 0.039 to form 4 images of the quasar. The 
images lie in the bulge of the galaxy, where the optical depth for microlensing is high and 
there is little dark matter. Furthermore, the low lens redshift leads to unusually short mi- 
crolensing variabi lity time scales (~ months). Microlensing fluctuations were first detected by 
Irwin et al.l (119891 ) and have now been monitored f or well over a decade, principally by the Qp - 
tical Gravitational Lensing Experiment (OGLE) (jWozniak et al.ll2000l : lUdalski et al.l 120061 ). 
Moreover, due to the symmetry of the system, and the proximity of the lens galaxy, time de- 



lay c orrections are negligible (< 1 day, iDai et al.ll2003l : IVakulik et al.ll2006l : iKoptelova et al. 



20061 ) which simplifies analysis. Because of these characterist ics, Q 2237-1-0305 has been the 



object of many studies, and monitored iii different bands (e.g..lWoznia. 



20031: 



Vakulik et al 



20061 : knguita et al.l[2008l: lEigenbrod et al. 



c et al 



2008a 



2000 



Dai et al. 



Mosquera et al 



2009 



O'Dowd et al.ll2010l ). Ou r X-ray studies ( Chen et al.l I2OIII . I2OI2I ). some of which were in- 
dependently analyzed by IZimmer et al.l (120111 ). not only found evidence for microlensing 
variability but also for a dependence upon the X-ray continuum energy. 



In this paper we present a complete analysis of the OGLE and lChen et al.l (l201ll . |2012| ) 
Chandra data to set constraints on the structure of the optical and X-ray emitting re- 
gions. We also compare our X-ray size estimates to those obtained for HE 435—1223 



(jBlackburne et al.ll2011bl ). QJ 0158-4325 (IMorgan et al.ll2012h . RXJ 1131-1231 (jPai et al. 



2OIOI ) and HE 1104—1805 (Blackburne et al. 2012), to explore possible correlations of the 
X-ray sizes with BH mass. A summary of the observations, a description of the lightcurve 
fitting method, and a discussion of the main results is presented in Section [2l In Section [3] we 
summarize the results and their implications. Throughout this work we assume Qm = 0.3, 
= 0.7, and Hq = 72 km sec"^ Mpc'^ 



2. Data Analysis and Discussion 
2.1. Observations 



Q 2237+0305 was obser ved in X-rays with AC IS (jGarmire et al.ll2003l ) on the Chandra 
X-ray Observatory ( Chandra AWeisskopf et al.ll2002l ). Combined with archival data, we have 
a total of 20 epochs between JD 2251 and JD 5529 taken during Cycles 1, 2, 7 and 11. The 
full band was defined as 0.4—8.0 keV (in the observers frame). These were then divided into 
soft (0.4—1.3 kev) and hard (1.3—8 keV) bands, where the division at 1.3 keV was chosen 
to balance the counts between the bands. A d etailed descriptiori of th e X-ray observations 
and data reduction techniques can be found in IChen et al.l (|201ll . |2012| ). 
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The optical data consist of the V-band observations from OGLE. The hghtcurve for 
our analysis goes from JD 2085 (June 2001) to JD 4602 (May 2008), and consists of 247 
ep ochs. Details of the OGLE data reduction and photometric techniques can be found 
(|2006[ ). We include the estimates of additional systematic errors in the 



m 



Udalski et al 



photometry of 0.02, 0.03 , 0.04 an d 0.05 mag for images A, B, C and D respectively, following 
Poindexter fc Kochanekl J2010al lbll Figures [1] and [2] show the flux ratios for the B, C, and D 



images relative to A. In Fig. |3]we have zoomed in 2006 (Cycle 7), where the X-rays show a 
significant variation over a short period of time, possibly due to a caustic crossing event. 



2.2. Magnification Patterns and Light curve fitting 



To fit the m icrol ensing lightcurves we applied th e Bayesian Monte Carlo Method of 
Kochanekl (12004r) and iPoindexter fc Kochanekl (j2010al ). closely following the procedures in 
Blackburne et al.l (j2011bl ). The basic idea is as follows. We built magnification patterns 
and convolved them with a source model to produce simulated light curves. We compared 
these light curves to the data and then use a Bayesian analysis to determine the parameters 
and their uncertainties. When fitting multiple-wavelength data sets, we first fit the best 
sampled lightcurve (OGLE V-band in our case), and use the same trials to sequentially fit 
the other bands (hard and soft X-rays here) assuming a new source size for each band. We 
next describe these steps in greater detail. 

For the models of the optical data we generated magnification patterns for each image 
of 4096x4096 pixels, with an outer scale of 20 x 20 Re, where for Q 2237+0305 the Ein- 
stein radius is Re = 1.8 x 1O^^(M/M0) ^'^cm. T his gives a resolution of 0.005 i?£;/pixel 
= 9.0 X 10^^(M/Mq)^/^ cm/pixel. For comparison, iMorgan et al.l ( 20101 ) estimated a single 
epoch CIV BH rn ass of lO^'^^M© with some caveats about the quality of the spectrum, and 
Assef et al.l ( 12011! ) estimated masses of lO^'^M© and 10^ '^ Mq for Ha and H/3 respectively, so 



the gravitational radius of the BH is roughly = GMbh/c ~ 1.5 x 10 (M^/z/IO^Mq) cm, 
and in this paper we eventually find an optical disk scale length of ^ 1.5 x 10^^ cm. We used 



the sam e surface density, k, and shear, 7, values as were used by iPoindexter fc Kochanek 
(l2010a( ). and fixed the mean microlensing mass to (M) = O.3M0, since this is roughly cor- 
rect for older stellar populations, and any small variation in mass would only produce a 
constant shift oc (M)^/^ in the marginalized distribution of the projected size. Since for this 
system the surface density is dominated by the stars rather than by dark matter, we as- 
sumed that all the surface mass densi ty is in stars, with k^/k = 1. In their earlier models of 
Q 2 237+0305 varying these quantities [Poindexter fc Kochanekl ( l2010al ) found (M) ~ O.3M0, 
and Kochanekl ( 120041 ) found ~ 1, consistent with these expectations. 
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The patterns change with time due to the motions of the Earth, the lens and its stars, 
and the source. For the Earth we combine the hehocentric CMB dipole and Earth's or- 
bit projected or ito the lens plane. We used th e source and lens galaxy peculiar velocities 
estimated from iMosquera &: Kochanekl (l2011bl ). We assigned the stars a random velocity 
dispersion of 170 km s~^ and velocity rotations o f 55 km s~^ for images A and B, and of 20 
km s~^ for images C and D based on lTrott et al.l (120101 ). 



Finally, the magnification pattern for each image and epoch is convolved with a simple 
thin disk model ( IShakura fc Sunyaevlll973l ) without including an inner edge. We adopted a 



simple profile since micro lensing is much raore sensitive to the source s ize than to the shape 
of the brightness profile (IMortonson et al.ll2005l : ICongdon et al.l 120071 ). and we ignored the 



inner edge because it is generally small compared to the measured (optical) disk sizes. The 
optical disk scale length, Rx, is defined by the radius at which ksTeffij-) = hc/X, where A 
is the observed wavelength. The X-rays, however, are better characterized by the half-light 
radius, R1/2 = 2.44i?A, and for simpli city we used the s ame b right ness profile since ther e 
seems no reason to change it based on IMortonson et al.l (120051 ) and ICongdon et al.l (120071 ) . 
While microlensing can constrain disk s tructure, it is generally achieved by rneasur ements 
of size as a function of wavelength (e.g., iMorgan et al.ll2012l : iBlackburne et al.ll2012l ) rather 
than distinguishing between detailed structures at fixed wavelength. 

Since the OGLE V-band light curve is the best sampled, consisting of 247 epochs from JD 
2085 to JD 4602, we fit it first. We traced 10^ random trajectories for each source structure, 
parametrizing the source models by choosing from projected areas covering logiQ{siTea./ cm^) 
= 29.4 to 34.6 in steps of 0.2, at each of 10 inclinations, i, with cos i values going from 
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cted substructure (e.g.. 


Falco et al. 


1999; 


Eieen 


3rod et al. 


2008b; 




et al.l 


2009), undete 


Mao & Schneider 


1998 




Metcalf & Madau 


2001 




Kochanek & Dalai 


2004; 
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2OI0I). uncertainties in the 


macro models, and 



any other contamination from the lens or host galaxy, the time-constant relative magnifica- 
tions of the images were allowed to vary with a Gaussian prior set to a 1.0 mag dispersion. 

We saved all the parameters for trials with x^/^dof < 5. For each of these trial pa- 
rameters we recompute the lightcurves for the next band, hard X-rays, but again varying 
the source size (reduced to a log^g (area/cm^) range of 28.8 to 34.2), inclination angle, and 
relative magnification offsets. Since the X-ray emission is more compact, we doubled the 
resolution of the patterns to 8192 x 8 192 pixels while rnaintairi ing the outer scale of 20Re 
following the procedures developed by IBlackburne et al.l (l2011bl ). This gives a pixel scale of 
2.4 X lO^^cm/pixel, which is small enough to resolve the average sizes found in the runs. We 
choose to fit the hard band next because its variabihty amplitude, especially during 2006, 
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appears to be larger than observed in the soft band. This will reduce the computational time 
by reducing the number of solutions found at this stage, and therefore the number of initial 
conditions for the next band. In later runs of this band sequence we can be less "liberal" 
in keeping solutions. We keep hard X-ray fits with l^dof < 4 and then fit the soft X-ray 
band saving results with x^ /Ndof < 3. Dropped solutions would contribute to our Bayesian 
probability integrals as exp[— (x^ — Xmj„)/2], so the dropped cases should be exponentially 
unimportant to the final results. 

Once the initial set of trials has been "filtered" through the three bands, we use a 
Bayesian analysis to determine the parameters and their uncertainties. Figure H] shows the 
probability distribution for the scale radii obtained with a logarithmic prior on the sizes. 
We found that the X-ray emitting regions are significantly smaller than the optical, and 
found no significant difference between soft and hard bands, although the uncertainties are 
large. Between the different bands we found half-light radius ratios of \og{Ri/2,softl Rv) = 
-0.651°:^^, \og{Ri,2Mrd/ Rv) = -0.5li°:^^, and \og{Ri,2,soft/ Ri/2Mrd) = -0.14+°:?°. We 
also ran a sequence considering only the OGLE V-band and the full X-ray band, since the 
X-ray uncertainties are smaller when considering the complete energy range. We saved trials 
with full X-ray band fits with l^dof < 5, and we found \og{Rxjuul Rv) = — 0.52+q;54 (Fig. 
|5]). Thus the X-ray half-light radius is on scales of ~ xO^'^+°'^ Vg, while the optical half-light 
radius is on scales of lO^'^+^'^r^ for Mbh = lO^'^M©. We also examined the distributions 
of the projected area ratios between the different wavelengths to see if the structure of the 
X-ray emitting region could be better defined (see Fig. These ratios have values of 
\og{Ax,fuu/Av) = -IMt'oil \og{Ax,soft/Av) = -1.21^?;^?, log(Ax,wMv) = -1-01^?:?^, 
and \og{Ax,soft/ Ax,hard) = — 0.21+q;9^ that are consistent with the derived half-light radius 
ratios. 

Despite finding large numbers of good fits to the V-band data alone, when we included 
the X-ray data we had difficulty obtaining large numbers of trials that were good fits to 
the datalj. Part of the problem is that it is probably expoi ientially harder to fit temporall y 



longer light curves using the Monte Carlo approach (see iPoindexter fc Kochanekl |2010aj) . 
But another problem could be that we do not have fully aligned optical and X-ray data 
because we have no extension of the OGLE light curves to the most recent X-ray epochs. 
We tested this hypothesis by fitting the hard band alone. The number of good solutions 
we found was again small, although no significant change in the size estimate was observed 
when compared to the multi-band run. This suggests that the problem lies in the length 
of the X-ray lightcurves. One of the effects of finding small numbers of solutions is that 



iBest ixyNdof )hard = 1.9, and best (xV^do/)so/t = 1.3. 
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the size at V-band is slightly shifted to larger values than found when fitting only the V- 
band data, because solutions with higher values dominate the weight after the multi-band 
run. Moreover, in the hard X-rays, many of the "good" solutions ignore the presence of the 
probable caustic crossing during 2006, whether only in the X-ray data or in the multi-band 
fits. 

To examine these questions, we considered all the parameters from the V-band only 
fits (the complete OGLE light curve), but then separately analyzed different portions of 
the X-ray light curves. In one sequence we only fit the caustic crossing event occurring 
during 2006 (from JD 3745 to JD 4115), and in a second sequence we fit the points that 
characterize the long time scale changes between the Chandra cycles (JD 1793 and 2251 
in 2000-2001, JD 3745 in 2006, and all the points in 2010, JD 5197 to 5527). Although 
analyzing them separately is less powerful in a statistical sense, we will not lose either the 
information encoded in the caustic crossing, or in the variability observed from cycle to cycle. 
Moreover higher quality solutions (i.e., with lower x^) will be obtained for each multi-band 
run because we have decoupled the "global" statistics on long time scales from the "local" 
statistics of a particular caustic crossing event. We otherwise followed the same analysis 
procedures. 

The probability distributions for the scale radii for the independent sequences are shown 
in Figure [3 From the caustic crossing fitting results (2006 data. Figure [71 top panel) good 
constraints can be set on the hard X-ray emitting regioij^ (log(i?i/2, hard/cm) = 15.411q;27), 
and as expected, due to it flatness, the soft band data does not provide any information 
about the size of the emission region in this energy range. The sparse long term X-ray 
data (Figure [TJ middle panel), constrains both the soft and hard X-ray emitting regionqfl, 
although the uncertainties are large (log(i?i/2, hard/cm) = 15.75lo;48 and log(i?i/2,so/t/cm) = 
15.70^035) • We can combine these results to improve our measurements by multiplying 
the probability distributions of the independent sequences (Fig. [71 bottom panel). For 
logarithmic size priors we find log(i?i/2,wd/cm) = 15.46+o;29, log(i?i/2,so/t/cm) = 15.76l:o:34, 
and log(i?i/2,v'/cm) = 16.41 ± 0.18. These results are consistent with the analysis of the full 
lightcurves, but we view them as more reliable given that we found a considerably higher 
number of good solutions for each sequence. 



^For the caustic crossing sequence the best {x^ /Ndof)hard ~ 0.9. 

•^For the sparse long term sequence the best {y^ /Ndof)soft — 1.02, and the best {x^ / Ndoj)hard = 1-3. 
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3. Summary and Conclusions 



M otiva ted by tlie recent detection of energy- dependent mic rolensing by IChen et a 



fl201ll . I2OI2I) in Q 2237+0305, we analyzed the OGLE V-band and lChen et al.l (120 111 120121 ) 



Chandra data to study the structure of the optical and X-ray emitting regions of this 
lensed quasar. We found that X-ray emission regions are significantly smaller than the 
optical, and considering two energy bands in the X-rays we estimated the half-light radii 
to be log(i?i/2,w/cm) = 15.46^0:29, log{Rx,soft/cm) = 15.76lo;34, and log(i?i/2,y/cm) = 
16.41 ± 0.18. These results also suggest a structure for the X-ray emission region in which 
the hard component comes from a more compact region, although the statistical signifi- 
cance of this result is low. Since the geometry and extension of the disk depend upon the 
physics that triggers the X-ray emission, generally attributed to inverse Compton scattering 
of soft U V photons from the accret ion disk by hot electrons in a corona surrounding the 



disk (e.g., [Reynolds fc NowakI 120031) . the smaller X-ray extent favors geor aetrical configu- 



ratio ns like the light bending ( iFabian et al.l 120051 ) or aborted jets models (IGhisellini et al. 
20041) and also agrees with th e general relativistic MHD models of iHirose et al.l (120041 ) and 
Machida fc Matsumotol J^OoJ. 



We now have enough measurements to consider correlations between the X-ray size 
and M bh- Fig. |H] shows all the available optical and X -ray sizes following Morgan et al. 



Fohlmeister et al.l mm . lOai et al.l fl2010). and 
results (open symbols) are from lDai et al.l ( I2OIOI ). 



Blackburne et al 



Morgan et al 



teOllbl. I2OI2 ), and t he B H mass estini a tes ar e from lPeng et al 



(I2OI0I). Here the o ptica l results (fil l ed squ ares) are from iMorgan et al.l (l2006l. I2OIOI |2012|) 



(l2011bl. l2012h. :he X-ray 



( 2012 ). andlBlackburne et al 



f l2006l ). Morgan et al.l f l2010h 



Greene et al.l (|2010[ ) , and lAssef et al.l (120 111 ) . Although X-ray size measurements exist for 
only a small number of systems, a first inspection suggests that the size of the X-ray emission 
region is roughly proportional to the mass of the central BH. Since the optical s izes appear 



to scale as the RopUcai oc predicted by thin disk theory ( iMorgan et al.ll2010l ). while the 



X-ray sizes appear to scale as Rx oc Mbh, the differences bet ween microlensi ng at X-ray and 
optical should be largest at low mass (e.g., RXJ 1131— 1231. iDai et al.ll2010l ) and modest at 
high mass (e.g., Q 2237+0305). However, the prevalence of upper limits on the X-ray sizes, 
as well as the large error bars, does not allow us to set strong constraints on the correlation. 
Because the X-ray emission is so compact, robust lower limits on sizes depend on having 
better sampled X-ray light curves than are typical of the existing data, although we are 
addressing this problem in a new set of Chandra observations. With additional and better 
measurements, we can not only better constrain the correlation with Mbh but also examine 
whether the spatial structure of the corona is correlated with the X-ray spectral index, and 
explore the origin o f the reflected X-ray components, particularly the Fe Ka lines, as we 
started to explore in IChartas et al.l (120 12l ). 



-9- 



Acknowledgments: This research was supported by NASA/SAO grants GO6-7093X, GOO- 
11121A/B/C, and G01-12139A/B/C, and NSF grants AST-0708082 and AST-1009756. Fur- 
ther support for this work was provided by the National Aeronautics and Space Adminis- 
tration through Chandra Award Number 11121 issued by the Chandra X-ray Observatory 
Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf 
of the National Aeronautics Space Administration under contract NAS8-03060. 



- 10 - 



REFERENCES 

Abajas, C, Mediavilla, E., Muoz, J. A., Popovic, L. C, & Oscoz, A. 2002, ApJ, 576, 640 

Agol, E., Gogarten, S. M., Gorjian, V., & Kimball, A., 2009 ApJ, 697, 1010 

Anguita, T., Schmidt, R. W., Turner, E. L., Wambsganss, J., Webster, R. L., Loomis, K. 
A., Long, D., & McMillan, R. 2008, A&A, 480, 327 

Assef, R. J., et al. 2011, 742, 93 

Bate, N. F., Floyd, D. J. E., Webster, R. L., & Wyithe, J. S. B. 2011, ApJ, in press 
flarXiv: 1103. 03531) 

Bate, N. F., Floyd, D. J. E., Webster, R. L., & Wyithe, J. S. B. 2008, MNRAS, 391, 1955 

Blackburne, J. A., Kochanek, C. S., Chen, B., Dai, X., & Chartas, G. 2012, in prep 

Blackburne, J. A., Kochanek, C. S., Chen, B., Dai, X., & Chartas, G. 2011, ApJ, submitted 
f!arXiv:1112.0027p 

Blackburne, J. A., Pooley, D., Rappaport, S., & Schechter, P. L. 2011, ApJ, 729, 34 
Blackburne, J. A., & Kochanek, C. S. 2010, ApJ, 718, 1079 

Blaes, O. M. 2004, in Accretion Discs, Jets and High Energy Phenomena in Astrophysics, 
Vol 78, ed. Vassily Beskin et al. (Berlin: Springer), 137 

Blaes, O. 2007, ASPC, 373, 75, "The Central Engine of Active Galactic Nuclei", ed. L. C. 
Ho and J.-M. Wang, (San Francisco: ASP) 

Chartas, G., Kochanek, C. S., Dai, X., Moore, D., Mosquera, A. M., & Blackburne, J. A. 
2012, ApJ, 757,137 

Chartas, G., Kochanek, C. S., Dai, X., Poindexter, S., & Garmire, G. 2009, ApJ, 693, 174 

Chartas, G., Agol, E., Eracleous, M., et al. 2002, ApJ, 568, 509 

Chen, B., Dai, X., Kochanek, C. S., et al. 2012, ApJ, 755, 24 

Chen, B., Dai, X., Kochanek, C. S., et al. 2011, ApJ, 740, 34 

Congdon, A. B., Keeton, C. R., & Osmer, S. J. 2007, MNRAS, 376, 263 

Dai, X., Chartas, G., Agol, E., Bautz, M. W., & Garmire, G. P. 2003, ApJ, 589,100 



Dai, X., & Kochanek, C. S. 2009, ApJ, 692, 677 

Dai, X., Kochanek, C. S., Chartas, G., Kozlowski, S., Morgan, C. W., Garmire, G., & Agol, 
E. 2010, ApJ, 709, 278 

Eigenbrod, A., Courbin, P., Dye, S., Meylan, G., Sluse, D., Vuissoz, C., & Magain, P. 2006a, 
A&A, 451, 747 

Eigenbrod, A., Courbin, P., Meylan, G., Vuissoz, C., & Magain, P. 2006b, A&A, 451, 759 
Eigenbrod, A., Courbin, P., & Meylan, G. 2007, A&A, 465, 51 

Eigenbrod, A., Courbin, P., Meylan, G., Agol, E., Anguita, T., Schmidt, R. W., & Wambs- 
ganss, J. 2008, A&A, 490, 933 

Eigenbrod, A., Courbin, P., Sluse, D., Meylan, G., & Agol, E. 2008b,A&A, 480, 647 

Elvis, M., 2000, ApJ, 545, 63 

Pabian, A. C, Miniutti, G., Iwasawa, K., & Ross, R. R. 2005, MNRAS, 361,795 
Palco, E. E., et al. 1999, ApJ, 523, 617 

Ployd, D. J. E., Bate, N. P., & Webster, R. L. 2009, MNRAS, 398, 233 
Pohlmeister, J., Kochanek, C. S., Palco, E. E., et al. 2008, ApJ, 676, 761 
Garmire, G. P., Bautz, M. W., Nousek, J. A., & Ricker, G. R. 2003, SPIE, 4851 
Greene, J. E., Peng, C. Y., & Ludwig, R. R. 2010, ApJ, 709, 937 
Ghisellini, G., Haardt, P., & Matt, G. 2004,A&A, 413, 535 

Guerras, E., Mediavilla, E., Jimenez- Vicente, J., Kochanek, C. S., Muoz, J. A., Palco, E.,& 
Motta, V. 2012, ApJ, submitted flarXiv:1207.2042p 

Hinshaw, G. et al. 2009, ApJS, 180, 225 

Hirose, S, Krolik, J. H., De Villiers, J. P., & Hawley, J. P. 2004, ApJ, 606, 1083 

Huchra, J., Gorenstein, M., Kent, S., Shapiro, I., Smith, G., Horine, E., & Perley, R. 1985, 
AJ, 90, 691 

Irwin, M. J., Webster, R. L., Hewett, P. C, Corrigan, R. T., & Jedrzejewski, R. I. 1989, AJ, 
98, 1989 



- 12 - 



Kelly, B. C, Bechtold, J., & Siemiginowska, A. 2009, ApJ, 698, 895 

Kochanek, C. S. 2004, ApJ, 605, 58 

Kochanek, C. S., & Dalai, N. 2004, ApJ, 610, 69 

Kochanek, C. S. 2006, in Saas-Fee Advanced Course 33, Gravitational Lensing: Strong, 
Weak and Micro, ed. G. Meylan, P. Jetzer, & P. North (Berlin: Springer), 453 

Koptelova, E. A., Oknyanskij, V. L., & Shimanovskaya, E. V. 2006, A&A, 452, 37 

Kozlowski, S, Kochanek, C. S., Udalski, A. et al. 2010, ApJ, 708, 927 

Lewis, G. F., & Belle, K. E. 1998, MNRAS, 297, 69 

MacLeod, C. L., Ivezic, Z., Kochanek, C. S., et al. 2010, ApJ, 721, 1014 

Mao, S., & Schneider, P. 1998, MNRAS, 295, 587 

Machida, M. & Matsumoto, R. 2003, ApJ, 585, 429 

Mediavilla, E., et al. 2009, ApJ, 706, 1451 

Mediavilla, E., et al. 2011, ApJ, 730, 16 

Metcalf, R. B., & Madau, P. 2001, ApJ, 563, 969 

Miller, L., Turner, T. J., Reeves, J. N. 2008, A&A, 483, 437 

Morgan, C.W., Kochanek, C. S., Morgan, N. D., & Falco, E. E. 2006, ApJ, 647, 874 

Morgan, C. W., Kochanek, C. S., Dai, X., Morgan, N. D., & Falco, E. E. 2008, ApJ, 689, 
755 

Morgan, C. W., Kochanek, C. S., Morgan, N. D., & Falco, E. E. 2010, ApJ, 712, 1129 

Morgan, C. W., Hainline, L. J., Chen, B., et al. 2012, ApJ, 756,52 

Mortonson, M. J., Schechter, Paul L., & Wambsganss, J. 2005, ApJ, 628, 594 

Mosquera, A. M., Muiioz, J. A., & Mediavilla, E. 2009, ApJ, 691, 1292 

Mosquera, A. M., Mufioz, J. A., Mediavilla, E., & Kochanek C. S. 2011, ApJ, 728, 145 

Mosquera, A. M., & Kochanek, C. S. 2011b, ApJ, 738, 96 



- 13 - 



Mufioz, J. A., Mediavilla, E., Kochanek, C. S., Falco, E. E., Mosquera, A. M. 2011, ApJ, 
742, 67 



Motta, v., Mediavilla, E., Falco, E.,& Munoz, J. A. 2012, ApJ, in press flarXiv:1206.1582p 
O'Dowd, M., Bate, N. P., Webster, R. L., Waytli, R., & Labrie, K. 2010 f!arXiv:1012.3480p 
Ofek, E., O., Rix, H., W., & Maoz, D. 2003, MNRAS, 343, 639 

Peng, C. Y., Impey, C. D., Rix, H. W., Kochanek, C. S., Keeton, C. R., Palco, E. E., Lehar, 
J., & McLeod, B. A . 2006, ApJ, 649, 616 

Peterson, B. M. 1993, PASP, 105, 247 

Poindexter, S., Morgan, N., & Kochanek, C. S. 2008, ApJ, 673, 34 
Poindexter, S., & Kochanek, C. S. 2010a, ApJ, 712, 658 
Poindexter, S., & Kochanek, C. S. 2010b, ApJ, 712, 668 

Pooley, D., Rappaport, S., Blackburne, J., Schechter, P. L., Schwab, J., & Wambsganss, J. 
2009, ApJ, 697, 1892 

Pooley, D., Blackburne, J. A., Rappaport, S., & Schechter, P. L. 2007, ApJ, 661, 19 
Reynolds, C. S., & Nowak, M. A. 2003, Phys. Rep., 377, 389 
Rozanska, A., & Czerny, B. 2000, A&A, 360, 1170 
Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Siemiginowska, A., Bechtold, J., Aldcroft, T. L., McLeod, K. K., & Keeton, C. R. 1998, 
ApJ, 503, 118 

Sluse, D., Surdej, J., Claeskens, J. P., Hutsemekers, D., Jean, C, Courbin, P., Nakos, T., 
Billeres, M., & Khmil, S. V. 2003, A&A, 406, 43 

Sluse, D., Claeskens, J.-F., Hutsemekers, D., & Surdej, J. 2007, A&A, 468, 885 

Sluse, D., Schmidt, R., Courbin, P., Hutsemekers, D., Meylan, C, Eigenbrod, A., Anguita, 
T., Agol, E., & Wambsganss, J. 2011, A&A, 528, 100 

Sobolewska, M. A., Siemiginowska, A., & Zycki, P. T. 2004a, ApJ, 617, 102 

Sobolewska, M. A., Siemiginowska, A., & Zycki, P. T. 2004b, ApJ, 608, 802 



- 14 - 



Trott, C. M., Treu, T., Koopmans, L. V. E., & Webster, R. L. 2010, MNRAS, 401, 1540 

Udalski, A., et al. 2006, Acta Astron., 56, 296 

Vakulik, V., Schild, R., Dudinov, V., et al. 2006, A&A, 447, 905 

Vegetti, S., Koopmans, L. V. E., Bolton, A., Treu, T., & Gavazzi, R. 2010, MNRAS, 408, 
1969 

Wambsganss, J. 2006, in Saas-Fee Advanced Course 33, Gravitational Lensing: Strong, Weak 
and Micro, ed. G. Meylan, P. Jetzer, & P. North (Berlin: Springer), 453 

Weisskopf, M. C., Brinkman, B., Canizares, C., et al. 2002, PASP, 114, 1 

Wozniak, P. R., Alard, C., Udalski, A., et al. 2000, ApJ, 529, 88 

Zimmer, P., Schmidt, R. W., & Wambsganss, J. 2011, MNRAS, 413, 1099 



This preprint was prepared with the A AS macros v5.2. 



- 15 - 



2002 2004 2006 2008 2010 




2000 3000 4000 5000 

HJD - 2450000 



Fig. 1. — Optical (V-band, filled squares) and X-ray (0.4 — 8.0 keV, open squares) microlens- 
ing light curves for Q 2237+0305. 
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Fig. 2. — Optical and X-ray microlensing light curves for Q 2237+0305. Filled squares 
represent V-band, and open squares and triangles soft (0.4—1.3 keV) and hard X-ray (1.3—8.0 
keV) respectively. 
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Fig. 3. — Closer view of the Chandra light curves during Cycle 7. The OGLE V-band data 
overlapping this period of time are also shown for comparison. Filled squares represent 
V-band, and open triangles and squares hard and soft X-ray bands respectively. 
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Fig. 4. — Probability distributions for the half-light radii of the V-band (solid), soft X-ray 
(dash-dotted), and hard X-ray (dashed) emission. The gray area correspond s to predicted 



value s of Vg, where the lower and upp er limits are for M rh = 0.9x lO^M© (CIV. lMorgan et al. 



20101) and Mbh = 2.4 x lO^M© (H/3, lAssef et al.ll201lh respectively. 
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Fig. 5. — Probability distributions for the half-light radii of the V-band (solid) and full X-ray 
(dashed) emission. The gray area corresponds to predicted values of r^, where the lower and 
uppe r limits are for M rh = 0.9 x lO^M© fCIV. iMorean et al.lboiol ) and Mbh = 2.4 x IO^Mq 
(H/3. lAssef et al.ll201ll ) respectively. 
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Fig. 6. — Top panel: Probability distributions for the logarithm of the half-light radii hard 
X-ray area to optical area ratio (solid), along with similar distributions for soft X-ray to 
optical (dashed) and hard to soft X-ray (dash-dotted). Bottom panel: Similar distribution 
for the full X-ray area to optical area ratio. 
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Fig. 7. — Probability distributions of the half-light radii for the V-band (solid), soft X-ray 
(dash-dotted), and hard X-ray (dashed) emission considering different sections of the X-ray 
light curves. The top panel shows the results corresponding to the analysis of the 2006 
data, the middle panel to the sparse long term X-ray data, and the bottom panel shows the 
combined results. As in Fig. HI the gray area corresponds to predicted values of Vg. 
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Fig. 8. — Emission region half-light rad ii as a function o f black hole mass Mrh bas ed in 
part on results fromlchartas et all J2009h . lDai et all J201o[ ). iMorgan et all J201ol . 120121 ) . and 
Blackburne et al.l (l2011bl . |2012| ). Filled squares represent optical sizes, and open triangles 
and squares are for the hard and soft X-ray sizes, respectively. Small offsets between the 
soft and hard X-ray data were applied for clarity. Diamonds correspond to full band X- 
ray size estirn ates. Upper limits are indicated by arrows. The dashed line corresponds the 
Morgan et al. (12010 ) fit to the optical data. 



